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Abstract 

Low-energy Landau levels of AB-stacked zigzag graphene ribbons in the presence of a uniform 
perpendicular magnetic field (B) are investigated by the Peierls coupling tight-binding model. 
State energies and associated wave functions are dominated by the B-field strength and the kz- 
dependent interribbon interactions. The occupied valence bands are asymmetric to the unoccupied 
conduction bands about the Fermi level. Many doubly degenerate Landau levels and singlet curving 
magnetobands exist along kx and k^ directions, respectively. Such features are directly reflected 
in density of states, which exhibits a lot of asymmetric prominent peaks because of ID curving 
bands. The fc^-dependent interribbon interactions dramatically modify the magnetobands, such as 
the lift of double degeneracy, the change of state energies, and the production of two groups of 
curving magnetobands. They also change the characteristics of the wave functions and cause the 
redistribution of the charge carrier density. The fc^-dependent wave functions are further used to 
predict the selection rule of the optical transition. 
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I. INTRODUCTION 



Carbon-based materials, such as^raphite, carbon nanotubes,[l|] two-dimensional (2D) 



graphene and few-layer graphenes, 



7 ^7 



% Is] and graphene ribbons, have been widely 
studied both experimentally and theoretically. Graphene ribbons could be synthesized 
by using heat treatment, B, 7, 8] pulsed-laser deposition technique, ^ patterning epitaxi- 

mm 



ally grown graphenes, [lO|, lUj tailoring exfoliated graphenes 



12 



microscopy. 



13| by scanning tunneling 



14i \ or chemical vapor deposition. [15'] A one-dimensional (ID) graphene ribbon 



is obtained by cutting a 2D graphene, planar hexagonal lattices of carbon atom, along the 
longitudinal direction. As a result of the special geometric structure, graphene ribbons have 
motivated many interesting studies on magnetic properties, 16|, 



17 
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erties, Hj, m 



191 . 120| optical prop- 



2l| electronic excitations, 22j or electronic properties. [23| The objective of 



this work is to investigate the magnetoelectronic structures of Bernal graphene ribbons. 

ID zigzag and armchair graphene ribbons, sided with two parallel zigzag and armchair 
structures along the longitudinal direction, respectively, are intensively studied. Many theo- 
retical studies predict that a zigzag ribbon has peculiar edge states. Such edge states produce 
flatbands at low energy and give rise to a conspicuous peak in density of states. 2J| An arm- 
chair ribbon does not exhibit such states. Its electronic properties, such as the band gap, 
depend on the ribbon width. Furthermore, 2D multilayer zigzag (armchair) graphene ribbons 
are the stack of infinite identical ID zigzag (armchair) graphene ribbons along the stacking 
direction [figure 1] . This stacked system exhibits the anisotropic energy dispersions between 
in-ribbon-pane and the stacking direction for the intraribbon coupling is much stronger 
than the interribbon interaction. Moreover, the interribbon interactions also modify the 
in-ribbon-pane electronic properties, such as state energies, energy dispersions, band-edge 
states, and size of band gap. The study results show that the geometrical structures (the 
stacking types (AA or AB stacking), ribbon width, and ribbon edge structure) have a signifi- 



cant effect on the electronic and optical properties of multilayer graphene ribbons. 25 



26 



27| 



The optical measurement might serve a method to determine the geometrical structures of 
graphene ribbons. ^ 



When a monolayer graphene ribbon is submitted to a perpendicular magnetic field B, 
its magnetoelectronic properties, e.g., the magnetic bands and the related wave functions, 
are determined by the competition between the magnetic confinement and the quantum 
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confinement. 



19 



20 



291] When the spatial extent of Landau wave functions is smaller than 



ribbon width, i.e., the magnetic confinement predominates over the quantum confinement, 
the Landau states exist and the state energy follows E oc a/HB (n meaning the subband 
index). [2^ Such a siniple characteristic is absent from a bilayer Bernal graphene ribbon. Ac- 
cording to literature, [30, 



31 



331] bilayer graphene ribbons exhibit different B-dependent 



Landau levels as a result of the interribbon interaction. In the energy region \E\ < 50 meV, 
the Landau-level energies are linearly dependent on the magnetic field strength. In the 
higher energy region, they deviate from the B-dependence to the v^-dependence with the 
increase of the field strength. Based on the findings mentioned above, Bernal graphene rib- 
bons (AB-stacked graphene ribbons) are thus expected to have different magnetoelectronic 
properties due to the different geometrical structure and interribbon interactions. 

On the other hand, experimental measurements confirm that Landau levels of mono- 
layer graphene and epitaxial graphene obey the relation E oc a/]^]B, Landau levels of 



Dirac fermions. 
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35 ]_^ Interestingly, Landau levels of Dirac fermions in graphite are 



experimentally observed. 361] Recently, Far infrared magnetotransmission measurements on 
a thin graphite show the V^-dependent absorption lines at the H point of graphite. 371] The 
magneto-transmission of graphite at K point, where interlayer interactions are maximum 
along the HKH edge of Brillouin zone, exhibits S-dependent Landau levels. 371] Both the 
theoretical and experimental studies on the electronic and magneto-optical properties of 
graphitic systems motivate us to study the magnetoelectronic properties of Bernal graphene 
ribbons. The organization of the present paper is as follows. The analytic Hamiltonian ma- 
trix elements of the tight-binding method for magnetoelectronic properties of AB-stacked 
graphene ribbons are first derived in Sec. 2. Then, Sec. 3 investigates the effects of mag- 
netic fields and the interlayer interactions on the band dispersions. Landau plot, and wave 
functions, following by the prediction of the selection rule of the optical transition. Lastly, 
conclusions are drawn in Sec. 4. 



II. THEORY 

The AB-stacked zigzag graphene ribbon is chosen for the model study. It has hydrogen- 
terminated zigzag edges along the longitudinal (x) direction [Fig. 1(a)]. The ribbon width 
is defined by the number {Ny) of zigzag lines along the y axis. B and A atoms are different 



3 



sublattices of biparticles hexagonal lattice of a graphene ribbon. There are 2Ny carbon 
atoms in a primitive cell of a monolayer graphene ribbon. To obtain the 2D AB-stacked 
zigzag graphene ribbon, the same ID zigzag ribbons are piled in the ABAB sequence along 
the z axis with the stacking distance = 3.35 A. The A atoms have the corresponding 
atoms on the x — y plane directly above or below them, while the projections of the B atoms 
are located at the center of the hexagonal rings directly above or below them [Fig. 1(b)]. 

In the presence of a perpendicular magnetic field B = (0, 0, B), the electronic properties 
of AB-stacked zigzag ribbons are modeled by the Peierls coupling tight-binding model. The 
Hamiltonian is 

^ = Y.^l^,e^"^'-^ctc, + H.c.l (1) 

where the subindices i,j denote the summation over all sites in a primitive unit cell. cl{cj) 
is the creation (annihilation) operation, which generates (destroys) an electron at i (j) site. 



The details of the atom-atom interactions [381] 7,^- = 70, 7i,--- ,76 are described in Fig, 



1(b). e*^'^^''J is the Peierls phase shift due to the applied magnetic field. Oij is the line 
integral of vector potential A from i to j in a unit of the flux quantum $0 = ch/e. To keep 
the translation invariance along the x axis, within the Landau gauge, the vector potential 
A = (—By, 0, 0) is deliberatively adopted. In this gauge, AB-stacked zigzag ribbons have 
two ribbons or ANy carbon atoms in a primitive cell. The first Brillouin zone is a rectangle 
defined by — vr//^ ^ kx < tt/Ix, and —7t/2Iz < < 7r/2Jz. 

The tight-binding Bloch function is the linear combination of basis l-B^); l^m)> 

and l-B^)) which are the linear superposition of the 2pz orbitals located at A^, B^, A^, 
and B^ carbon atoms. The wave function is 

Ny Ny 

mkx,y,k,)) = J2 ^Ai^lAln) + bs^jBl) + + bBi\Bi), (2) 

m=l m=l 

where wave function coefficients a/ii , bai , ua'^ and 6^2 are the site amplitudes. 

The Hamiltonian representation is a 4Ny x ANy Hermitian matrix and could be regarded 

as 



H 



hi hi2 
h2i h2 



(3) 
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where each Hermitian block matrix /ii,/i2,^i2 or /i2i has 2Ny x 2Ny elements. The block 
matrix hi, Hamiltonian representation of the lower ribbon, is a tridiagonal matrix and its 
matrix elements hi{i,j) are: 

hii2m - 1, 2m - 1) = 76 + /3^^5 /2, 

hi{2m,2m) = p^-f 2/2, hi{2m - 1,2m) = hi{2m,2m-l) = 2-fo cos{V^bK/2 - n{m - [A^])$}, (4) 
hi{2m, 2m + 1) = hi{2m + 1, 2m) = 70, 

where m = 1, 2, ■ ■ ■ , A^^ and /? = 2cos{kzIz)- [N] = {Ny + l)/2 is used to locate the origin of 
coordinate in the center of the ribbon. 7r$, the Peierls phase shift, results from the magnetic 
flux passing half a hexagonal ring. The diagonal term 76 + /3^75 /2 (/3^72 /2), due to the 
interribbon interactions, can be treated as the equivalent site energy of A (B) atoms. The 
centers of the two ribbons in a unit cell do not coincide with each other [Fig. 1(a)]. This 
leads to a different tridiagonal matrix /i2, which is 

/i2(2m-l,2m-l) = 76 + P^^5 /2, 
h2{2m,2m) = (3^-f2/2, 

/i2(2m - 1, 2m) = /i2(2m, 2m - 1) = 270 cos{V^bk^/2 - 7r(m - [N] + i)<l>}, 
/i2(2m - 1, 2m + 2) = /i2(2m + 2, 2m - 1) = 70. 

The term 7r$/3 in Eq. (jSj) is caused by the difference between the centers of the two ribbons 
in a unit cell. 

The two off-diagonal block matrices satisfy a simple relation hu = /i2i- As a result of 
the interribbon interactions, the nonzero matrix elements huihj) of the off-diagonal block 
matrix follows: 

^/ii2(2m-l,2m-l) = /37i, 
hi2{2m,2m) = p'j-i, 

hi2{2m - 1, 2m) = 2/^74 cos{y/2,hk^/2 - 7r(m - [N] + |)<l>}, 
/ii2(2m, 2m - 1) = 2/^74 cos{^/?,hk^/2 - 7r(m - [N])^}, (6) 
hi2{2m, 2m + 1) = /574, 
/ii2(2m-l,2m + 2) = /374, 

/ii2(2m, 2m + 2) = 2/^73 cos{V3bk^/2 - n{m - [N] + i)$}. 



(5) 



The effect of magnetic fields is reflected in the interaction 27jCos(-\/3 &fcx/2 — n{m — [A^])$). 
All the nonzero elements of hi2 depend on (3 = 2cos{kzIz)- The interribbon interactions 
are equal to zero at k^Iz = 7i"/2, and they are maximum at k^Iz = 0. Eigenvalues and 
eigenvectors are obtained after the diagonalization of the Hamiltonian matrix. Eigenvalues 
are the energy dispersions E'^''"{kxlx, kzh), where c (f ) represents the unoccupied (occupied) 

; ; rpU • -J- I C,V 1C,V C,V 1C,V C,V 1C,V C,V 1C,V \ ■ ,1 7 

states. The eigenvector |a^i, o^i, ■ ■ ■ ^ ^ai,^ "bi,^ " ' ' "-Af^ "sf ' ' ' ^ ^a^^^ "b^^ " 
dependent envelope function \'^'^''"{y,kzlz)) along the y direction. 



III. MAGNETO-ELECTRONIC PROPERTIES 



The band structures of the 2D multilayer graphene ribbons exhibit a highly anisotropic 
structure. Figs. 2(a)-2(c) show the fc^-dependent low-energy bands for the AB-stacked 
zigzag graphene ribbon with width A^y=3000 (639 nm) at B = 0, B = 20 T, and B = 10 
T, respectively. The unit of energy is eV. The Fermi level is set to Ep = 0. At B = 0, 
the low energy bands are similar to those obtained from the LDA calculations. 391] These 
bands include parabolic bands and fiatbands [Fig. 2(a)]. The occupied valence bands are 
asymmetric to the unoccupied conduction bands about Ep = due to the energy-dependent 
interribbon interactions. The Ny = 3000 AB-stacked zigzag graphene ribbon is a semimetal 
because the valence bands lightly touch the conduction bands at Ep = 0. The degeneracy 
of the fiatbands at = 0, a special feature of a single zigzag ribbon due to the zigzag-edge 
boundaries, is lifted by the interribbon interactions. The partial fiatbands near Ep = 
correspond to the edge states mainly localized at the outmost zigzag positions. According 
to Figs. 2(b)-2(d), the magnetic field drastically modifies the energy bands, such as the 
alteration of the band feature, the shift of the subbands, and the production of the Landau 
levels. At B = 20 T, as shown in Fig. 2(b), the original parabolic bands might become the 
complete Landau levels. The occupied states E^ are asymmetric to the unoccupied states 



E'^ about Ep = 0. The Landau-level energies do not follow the simple relation E oc a/|?t-|B. 
The presence of B induces a longer and weak splitting energy subband near Ep, namely, the 
zero mode. There are more Landau subbands and a shorter range in band structure as the 
magnitude of the magnetic field decreases [Fig. 2(c)]. The zigzag boundaries make narrow 



graphene ribbons exhibit special electronic and magnetic properties. 



40 



Interestingly, the interribbon interactions induce the concave-upward and -downward 
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magnetobands along kz- The energy dispersions near Ep at B = 20 T with A^^ = 3000, B = 

10 T with Ny = 3000, and B = 20 T with A^'^^ = 6000 (1278 nm) along the fc^-axis are shown 
in Figs. 2(e)-2(h), respectively, to illustrate the effect of interribbon interactions. There 
are oscillating subbands, concave-upward and -downward bands, and nearly dispersionless 
bands nearby Ep [Figs. 2(e) and 2(f)]. The occupied states E^ are asymmetric to the 
unoccupied states E'^ about E = 0. The main feature of the energy dispersions remains 
unchanged for the variation of magnetic flux. However, the state energies are sensitive to 
the magnitude of the magnetic fleld [Fig. 2(g)]. The zero mode near Ep, including two 
oscillating subbands, has a dispersion of 0.04 eV at the low fleld limit. The upper subband 
of the zero mode, indexed by n = 0, shows the cosine energy bands E oc cos{kzIz) along 
the kz directions, while the feature of the lower subband, indexed by n' = 0, is distorted. 
A pseudogap A associated with the zero mode near Ep is observed. Thus, the size of the 
pseudogap A varies with kz. The upper limit of the pseudogap A, occurring at kzh = 7r/2, 
is equal to 70, the chemical difference between A atoms and B atoms. The concave-upward 
and -downward subbands, are classifled into two groups. One group, denoted as Group I, 
moves towards Ep while the other. Group II, moves away from Ep. The flrst three subbands 
of each group are shown in Fig. 2(f) as examples. The unoccupied states E^ and occupied 
states E^ of Group I are denoted by subband indices n and n', respectively. Those of Group 

11 are indexed, respectively, by n and n'. The subband crossings occur at different /c^'s. 
The band edge states, the local minimum and maximum, of the n, n, and n' subbands are 
located at kze = 0. 

The Landau plot, the Landau-level energies vs. the fleld strength, at different fc^'s can 
reveal the effect of the interribbon interactions. The Landau plot of the A^j, = 3000 AB- 
stacked zigzag ribbon at kziz = 7r/2 is shown in Figs. 3(a) and 3(b). The state energies of 
the zero mode are independent of the fleld strength. The energy spacing between the upper 
and lower subbands, the pseudogap A at kzh = 7r/2, is equal to 75. Except the zero mode. 
Landau-level energies follow the simple relation E^'"" oc a/B. The chief reason is that the 
two ribbons in a primitive unit cell are regarded as two independent ribbons because the 
interribbon interactions are turned off at kzIz = 7r/2. The magnetoelectronic properties are 
dominated by the magnetic conflnement rather than the quantum effect. Accordingly, the 
Landau-level energies, indexed by n, follow the relation E oc A^/|rj-|B. j2o| On the other hand, 
the interribbon interactions between the two ribbons are maximum at kzh = 0. The Landau 
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plots, as shown in Figs. 3(c) and 3(d), exhibit different behaviors. The state energies of the 
zero mode vary with the field strength. In the energy region \E\ < O.O270 ~ 50 meV, the 
state energies of the first and second unoccupied (occupied) Landau levels exhibit the linear- 
in-B dependence at i? < 50 T. The high energy Landau levels, indexed by n > 3, exhibit 
the linear-in-B dependence at 5 < 15 T. Then, they evolve from a linear B-dependence to 
a square-root B dependence with the increase in the field strength. Finally they exhibit the 
linear-in-A/B dependence at B > 20T. 

Density of state (DOS) of multilayer graphene ribbons is evaluated by D{uj) = Xlk/i 
|Vk-E''^(k; (f))\^i^^. As shown in Fig. 4(a), in the absence of a magnetic field, the featureless 
DOS of the Ny = 3000 graphene ribbon, asymmetric about u = 0, shows a small bump 
at = 0, which results from the tiny overlap between the occupied bands and unoccupied 
bands [Fig. 2(a)]. The A^^^ = 3000 graphene ribbon is a semimetal. The magnetic field 
modifies electronic states [Fig. 2] and the aspect of DOS [Fig. 4(b)]. DOS at B = 20 T 
is asymmetric about Ep. It chiefly exhibits three kinds of structures, including the delta- 
function-like peak, the compound peak near uj = 0, and the ID power-law divergences. The 
first kind of peaks are related to the Landau levels along [Fig. 2(b)] and dispersionless 
bands along [the red curves in Fig. 4(b)], the second kind to compound band along k^ 
and kz, and the third kind to concave-upward and -downward bands along kz [the red curves 
in Fig. 4(b)]. The sharp peaks at = —0.3 and —0.45 eV origin in the Landau levels along 
kx and dispersionless bands along kz. DOS of the concave-upward and -downward diverges 



in the forms 1/ a/cu — E{kze) and 1/ a/^(A!^^)~-^, respectively, where kze is the band edge 
of the subband. DOS exhibits more peaks in the decreases of the magnitude of magnetic 
fields, as shown in Fig. 4(c). In other words, the peak heights and positions are strongly 
dependent on the magnetic flux. Fig. 4(d) shows DOS of the A^j, = 6000 graphene ribbon 
at B = 20 T. Remarkably, density of states is independent of the ribbon width [Fig. 4(b) 
and Fig. 4(d)]. 

When the ribbon width is sufficiently large (A^ > 3000), the AB-stacked graphene ribbons 
at B = are expected to be almost identical to the bulk graphite in DOS except at very low 
frequency. The zigzag boundary induces partial fiat bands near Ep = and thus strongly 
affects the low-frequency physical properties. The two kinds of systems might differ from 
each other in the transport and magnetic properties. Nevertheless, the large graphene 
ribbons are useful in understanding the bulk graphite, e.g., magnetoelectronic properties [43| 
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and magneto-optical absorption spectra. The main reason is that the calculations of the 
magnetoelectronic structures are very complicated for the latter at B < 10 T (a very large 
Hamiltonian matrix). That is to say, the magnetoelectronic properties (state energies and 
associated wave functions) remain the same when the ribbon width is sufficiently wide. It 
is deduced to be the same with that of a 3D bulk graphite. 

The effect of interribbon interactions on the envelope functions deserves a closer exam- 
ination. Envelope functions of the Ny = 3000 zigzag ribbon subjected to B = 20 T at 
different fc^'s are shown in Figs. 5(a)-5(d). Each envelope function '^{y^kzlz) is separated 
into eight components, ij{Al), ij{Bl), ij{Al), ij{Bl), ij{Al), ^^(5^), ij{Al), and ^^(Ef). Al 
(-Bg) denotes A (B) atoms located at the odd (even) zigzag lines at the lower (upper) rib- 
bon plane. The subenvelope function ip{Al), for example, is |a^i, a^i, a^i, ■ ■ ■ ). Because 
^PiAl) = -ij{Al), ^{BD = ~^{Bl), ^{AD = -^{AD, and ^(5^) = -^(i?^), only four 
components ■ip{Al), ip{Bl), ip{Al), and 4'{B'^) are displayed. At kzh = 7r/2, the interribbon 
interactions are closed due to /? = 2cos{kzIz) = [Eqs. (|l])-([6])]. The two ribbons in a 
primitive unit cell are decoupled, i.e., there are two isolated graphene ribbons. Except the 
zero mode, the state energies are double degenerate. Their associated envelope functions, 
'^i{y,kzlz = 7r/2) and '^2{y,kzlz = 7r/2), are orthogonal to each other and different in 
the sign of some components. One of the two envelope functions is shown in Fig. 5(a). 



It is exactly described as 



MK)An-l{Bl),<Pn-l{Al),MBT. 



where 0„ is the harmonic 



oscillator, product of the Hermite polynomial if„ and the Gaussian function. 2^ Each 0„ 
has n nodes, where the subenvelope function changes the signs. The node number n of 
the subenvelope function iJjn{Al) is identical to the subband index n. Thus, we use the 
subband indices n = 1, 2, ■ ■ ■ and n' = 1, 2, ■ ■ • to specify the characteristic of envelope 
function [Fig. 5(a)]. (pniAl) and (pn^Bl) [(f)n-i{,Bl) and 0„,_i(v4g)] belong to the n \n — 1] 
mode. The envelope functions are the same as those of a single graphene ribbon. 20|] Due 
to the orthogonality, \E'„ and are different in the sign of some components. Besides, 
\E'„=o = 4>oiB'^) — 00 (-Be), the envelop function of the n = subband, is located in the B 
sites of the upper ribbon, (poi^l) — 4>o{^l), the envelop function of the n' = subband, 
is mainly controlled by the A atoms on the lower ribbon plane. The chemical difference 
76 changes the site energies of atoms A and, thus, causes the energy split and produces a 
pseudogap A = between \E'„=o and \E'„'=o at kzh = 2/n [Fig. 2(e)]. 

Away from kziz = 7i"/2, the switched-on interribbon interactions not only produce two 
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groups of subbands, but also change the aspects of envelop functions. At kzh — 0.4997r, 
the envelope functions related to n and n subbands concentrate on the lower and upper 
ribbon, respectively [Fig. 5(b)]. The subenvelope functions related to the n subband are 
(f)n{Al) and (j)n-i{Bl), where n = 1, 2, 3, ■ ■ ■ (n' = 1, 2, 3, ■ ■ ■). As for the n subbands, 
the associated subenvelope functions are (f)n{Al) and (j)n+i{B^), where n = 0, 1, 2, 3 - ■ •. 
Notably, n is different from n and it begins from n — {n' — 0). For example, the envelop 
function of n = 1 [n = 0] subband is *„=i = 0i(^^) +0o(5o) [*n=o = M^D+MB^)]- The 
effect caused by the interribbon interactions on the envelope functions can be interpreted 
by the first order perturbation theory. Around kziz — O.Stt, the matrix elements of hi2 
are proportional to (/t/2 — k^Iz) because of /9 = 2cos{kzIz) ~ 2(7r/2 — kzIz)- Thus, the 
weak interribbon interactions lift the degeneracy at kzIz = vr/2 and produces a two-level 
system, for example, n — 1 and n = subbands. The related envelope functions are 
'^{y,kzlz) = [^i{y,kzlz = 7r/2) ± ^2{y,kzlz = 7r/2)]/v^, the bonding or antibonding 
of the envelope functions at kzIz — 7r/2. As a result, the envelope functions of the n {n) 
subbands chiefly distribute on the lower (upper) ribbon with the different spatial symmetries. 
In addition, interribbon interactions have no influence on the envelop functions ^n=o and 

The feature of the wave function at kzh = 0.47r, as shown in Fig. 5(c), is dissimilar 
to that at kzIz — 0.4997r [Fig. 5(b)]. The charge carriers of the former distribute on two 
ribbon planes and each subenvelope function is of the same importance, whereas those of 
the latter concentrate only on one ribbon plane. The varying interribbon interactions causes 
the redistribution of charge density. The definition of subband indices n and n {n' and n') 
is based on the spatial symmetry of the subenvelope function ■?/'„(A^). The subenvelope 
functions are proportional to the harmonic oscillator 0„, for example, 7/'„=i(y4„) ~ (})i{A\). 
The envelope functions of the \n\ > 1 subband are [^^(^Iq), 0„_i(S^), 0„(74^), 0„+i(S^)] . 
The subband index n {n') is assigned to be n = 1, 2, 3, • • • {n' — 1, 2, 3, • • •). Notably, 
subenvelope functions ip{Al) and ■0(^o) belong to the n mode and ipiBl) [0(-B^)] to the n — 1 
[n-|- 1] mode. The spatial symmetry [n,n — l,n,n + l] is different from that [n,n— l,n — l,n] 
at kz = 0. On the other hand, the envelope functions of the n subband is the linear 
combination of components 0^(A^), (pn^i^Bl), 4>fi{Al) , and (f)n+i{B'^). The spatial symmetry 
[n,n — l,n,n + 1] of the n subband is the same as that of the n subband [n,n — l,n,n + l]. 
According to the spatial symmetry of ■0„(A^), the subband indices are assigned as ri = 



10 



0, 1, 2, • • • . Notably, n starts at zero. Due to the interactions between the n — ±6 and 
n — ±1 subbands, the shapes oi n — envelope functions is modified and its component 
ip{Bl) almost disappears. 

From kziz = OAtt to kziz = 0, the gradually enhancing interribbon interactions signifi- 
cantly modify the feature or the site amplitude of the n envelope functions, as shown in Fig. 
5(d). The maximum interribbon interactions, occurring at kzh — 0, dramatically change 
the main feature of the envelope function [heavy dots in Fig. 5(d)], which now can not 
simply be described by harmonic oscillator 0„. For comparison, those at k^Iz — O.Stt are 
shown by light dots in Fig. 5(d), whose main feature of the envelope functions is still similar 
to those at kziz = OAn. Thus, the envelope function is approximately linear superposition 
of components 0„(A^), 0„_i(-B^), 0„(A^), and 0„+i(i?^). The interribbon interactions dra- 
matically alter the site amplitude. ip!^\{Bl) and ip^!^i{Bg) dominate the envelope function. 
The comparison shows that the increase in the interribbon interactions not only change the 
main feature of the subenvelope functions but also enhance the site amplitude of B atoms. 
It should be noted that the interribbon interactions do not break the spatial symmetry of the 
envelope function, which is denoted as [ilJn{Al) , ipn-i{Bl) , ipn{Al) , ipn+iiB"^)] . The spatial 
symmetry of the envelope function has a critical effect on the optical selection rule. 

The dependence of magnetoelectronic properties of the n subbands on the wave vector 
kz is deliberated. For comparison purpose, the /c^-dependent state energies of the n — 1 and 
n — subbands and the associated envelop functions are shown in Fig. 6. At A;^ = 7r/2, two 
subbands are degenerate and magnetoelectronic properties mainly exhibit the characteristics 
of an isolated graphene because of the close of the interribbon interactions. As kz changes 
from 7r/2 to zero, the two degenerate subbands are lifted by the /c^-dependent interribbon in- 
teractions and form concave-upward and -downward bands along kz [Fig. 6(b)]. Meanwhile, 
the interribbon interactions have different influences on the shape of wave functions and the 
charge-carrier distribution [Fig. 6(a)]. The subenvelope functions oi n — subband [heavy 
dots in Fig. 6(a)] can be clearly described by the harmonic oscillator By contrast, the 
features of subenvelope functions related to n = 1 subband are seriously modified by the kz- 
dependent interribbon interactions as kziz 0. The /c^-dependent interribbon interactions 
cause the change of the site amplitudes. The charge carriers of the n = subbands tend to 
concentrate on A atoms, while those of n = 1 subbands are mainly located in B atoms [Fig. 
6(a)]. 
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The characteristics of wave functions are apphcable to determine the optical transition 
channels. The optical absorption spectra is A{lj) oc Yl^u Iuibz^ ^-^(^) ^ \^^l,v\'^■l 
where fi = n, n and u = n', n' are subband indices. D{u) is the joint density of states 



and Mfj,^^ is the velocity matrix element. 
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32| According to the result of Ref . [32] , the 



magnitude of velocity matrix is |M^,j,| = \Mf,^i,\intraribbon + |M^,i.|mter-rif>6on, where \Mf,^^\ is 
proportional to {'^^{y,kzlz)\'^u{y,kzlz)), the projection of the envelope function of initial 
state, l^^uiy, kziz)), onto that of final state, {"^^{y, kziz))- The subscript intraribbon [inter- 
ribbon) denotes the initial and final states located at the same (different) ribbons. The 
mtrari66on-projection, for example, is 

{'i/n{y,kzlz)\'^niy,kzlz)) intraribbon 

{MA')\ij,,iB')) + (8) 

The velocity matrix element M^,^ depends on the spatial symmetries of the initial 
and final states and, thus, determines the selection rule and the effective transition 



intraribbon^ 



channels. [19|, l20|, |42|] At kzh = v/2, only the intraribbon optical transition, |M^^j^| 
makes contribution to the velocity matrix element |M^^j,|. There is only one Group of sub- 
bands, i.e., n = n and z/ = n' . Therefore, the optical selection rule is predicted to be 
5n-n' = ±1, where = 0, 1, 2, 3 ■ ■ ■ and n' = 0, 1, 2, 3 ■ ■ ■ . The selection rule is the same 
as that of a monolayer graphene ribbon. 2o| 

Away from kzIz = 7i"/2, there are four possible optical transition channels n' n,n' n, 
n' —>■ n and n' —>■ n [Fig. 6(b)]. According to the spatial symmetry of envelope functions, 
the corresponding selection rules of the four transitions are predicted to be Sn-n' = ±1, 
6n-n' = ±1, Sn-n' = ±1, and 6n-n' = ±1- The change of the state energies caused by 
interribbon interactions might alter the peak positions of joint density of states, and the 
modification of the shape and site amplitude of envelope functions might give rise to the 
vivid variation of the magnitude of velocity matrix element Thus, the optical spectra 

of the AB-stacked ribbon are expected to be more complicated in peak number, peak position 
and peak height than those of a monolayer or bilayer ribbon. The work to calculate the 
optical spectra of the AB-stacked ribbon is on-doing. 
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IV. CONCLUSIONS 



The Peierls coupling tight-binding method is employed to investigate the anisotropic 
magnetoelectronic structures of multilayer graphene ribbons. They are strongly dependent 
on the B-field strength and the interribbon interactions. B can induce Landau levels along 
kx and change energy spacing. The interribbon interactions significantly affect state de- 
generacy, energy dispersions, Landau plot. They destroy the symmetry of magnetoband 
structures. They also produce concave and convex subbands along kz- The magnetoband 
structures are strongly anisotropic. There are many doubly degenerate Landau levels and 
singlet ID curving bands along k^ and k^ directions. Such features are directly reflected 
in density of states. DOS exhibits a lot of asymmetric prominent peaks, for example, the 
delta-function-like peak, the ID power-law divergence, and the compound peak. Moreover, 
the /c2-dependent interribbon interactions modify the shape of wave functions, alter the 
spatial symmetry, and arouse the change of charge-carrier distribution. The findings are 
further used to predict the optical selection rule. Most importantly, the large AB-stacked 
graphene ribbons are expected to be almost identical to the bulk graphite in low-energy 
magnetoelectronic properties. Our work is useful in understanding the bulk graphite, e.g., 
magnetoelectronic excitations and magneto-optical absorption spectra. The predicted mag- 
netoelectronic properties could be examined by transport and optical measurements. 
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Figure Captions 



FIG. 1. (a) The geometric structure of the Ny=3 AB-stacked zigzag graphene nanoribbons. (b) 
7o = 2.598el^ is the intralayer interaction. 71 = O.SGieV (75 = 0.036e\^) represents the 
interaction between two A atoms from two neighboring ribbons (two next-neighboring 
ribbons), and 73 = 0.319eV (72 =-0.014 eV) is for B atoms. 74 = 0.177eV corresponds 
to the interribbon interaction between A atoms and B atoms. 75 = — 0.026e\^ is the 
chemical-shift between A atoms and B atoms. The values of 7i's are the same as those 
of AB-stacked graphite, 



FIG. 2. The fc^;- dependent low-energy bands for the AB-stacked zigzag graphene ribbons with 
(a) iVj^=3000 at B = 0, (b) Ny=3000 at B = 20 T, (c) Ny=3000 at B = 10 T, 
and (d) Ny=6000 at B = 20 T. Figs, (e)-(h), same plot as Figs.(a)-(d), but for the 
/c^- dependence. 

FIG. 3. (a) Landau plot of the A^^ = 3000 ribbon at kziz = tt/2. (b) The Landau level energies 
of the Ny = 3000 ribbon at kj^ = tt/2 vs. VB. (c) Landau plot of the Ny = 3000 
ribbon at fc^ = 0. (d) The relation between Landau level energies of the Ny = 3000 
ribbon at /c^ = and \/B. 

FIG. 4. The low-energy DOS for the AB-stacked zigzag graphene ribbons with (a) Aj^=3000 
at B = 0, (b) Ny=3000 at B = 20 T, (c) ^"2^=3000 at B = 10 T, and (d) ^^^=6000 at 
B = 20 T. 

FIG. 5. The envelop functions for the AB-stacked zigzag graphene ribbons A^y=3000 subjected 
to B = 20 T at kzIz = 7r/2, 0.4997r, 0.47r and are shown, respectively, in (a)-(d). 

FIG. 6. (a) The energy dispersions of the n = 1 and n = subbands. (b) The variation of the 
envelop functions, related to the n = 1 and n = subbands, with k^Iz- 
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